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Abstract 

The time irreversible evolution of the two-nucleus spin systems interacting with a 
magnetic field is analyzed in the framework of the subdynamics theory. The spin 
systems are determined by the H(l) and C(13) nuclei. The investigation is based 
on the complex spectral representation of the Liouvillian which gives a rigorous 
description of irreversibility. The evolution of the density matrix is analyzed for the 
short and long time region. Quantum Zeno effect is discussed. 
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1 Introduction 



The irreversible nature of the world impels us to develop theory by adequate 
means, in accordance with this fact. However, the exact description of irre- 
versibility on the basis of fundamental classical and quantum theories is impos- 
sible, since they are reversible in the time. As a rule, the problem is solved on 
the basis of different physical approximations such as the Pauli master equa- 
tion [T], the Lindblad axiomatic formalism [2]. In these cases coarse-graining 
approximations, the distinction between an open system and its environment 
are used. The environment is assumed to be in thermodynamic equlibrium. 
Also, the problem is solved on the basis of the Schroedinger equation with a 
dissipative term [3], on the basis of the mixed quantun-classical approaches [I]. 
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On the whole, the approaches mentioned above, having a nature of the ap- 
proximations, contain a subjective, anthropomorphic element and cannot be 
completely satisfactory. Thus, the problem of describing the irreversible world 
on the basis of the reversible equations of classical and quantum physics i.e., 
the determination of the objective laws of nature arises. 
It is obvious that fundamental solution to this problem can be achieved when 
it is based on the equations which originally have irreversible nature. In this 
case, one should speak about the formulation of the dynamics, which makes 
it possible to include irreversibility in a natural way. 

For many years, the Brussels-Austin group, that was headed by I. Prigogine, 
has been developing an alternative approach with the aim to reveal the role 
of the irreversible processes, as well as their fundamental role in nature. A 
conclusion has been made that the solution to this task may be found outside 
the Hilbert space. In this approach, the eigenvalues problem is reduced to the 
complex eigenvalues problem [5], [5] 



H | ip a >= Z a | ip a >, < ip a | H =< ip a | Z, 



a- 



where the right eigenstate | <p a > and the left eigenstate < ip a | must be 
distinguished and the eigenvalue Z a is a complex number 

Z a = E a - ij a . (2) 

Here E a is a renormalized energy and 7 Q is a real positive value. The index a 
corresponds to the certain quantum state. In this case, the eigenstates | ip a >, 
< (p a | have a broken time symmetry. The eigenstate | <p a > is associated with 
the unstable state which vanishes for t — > +00, while | (p a > corresponds to 
the state vanishing for t — > —00 

I <p a (t) >= exp(-iE a t - j a t) I tp a (0) >, (3) 

I <p a {t) >= exp{-iE a t + la t) I £ Q (0) > . (4) 

The eigenstates | ip a >, < (p a \ are named the Gamow vectors and are in detail 
discussed in the literature [7] - [TO] . 

Further development of the Brussels-Austin group approach was obtained on 
the basis of the Liouville space extension of quantum mechanics [5] , [11] , in the 
framework of the complex spectral representation of Liouville-von Neumann 
operator L (Liouvillian) and subdynamics theory. In this case, Liouvillian 
has complex eigenvalues that break time-symmetry without introducing any 
approximations. The basic ideas of the approach, in detail, are examined in the 
works [5], [6], [11] - |14] . Here, I briefly mention some specific features. Let me 
examine the equation which determines the density operator p - Liouville-von 
Neumann equation 

i^ = M<). (5) 
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Liouvillian L has the form 

L = H x 1 - 1 x H, (6) 

where symbol " x" denotes the operation (A x B)p=ApB, and H is the total 
Hamiltonian, with H = H + Hj being the sum of the free part H and the 
interaction part Hj. Similarly, the operator L can be written down as the sum 
of the free part Lq, that depends on the Hamiltonian Hq and the interaction 
part Lj, which depends on Hi : L = L + Lj. Usually, for the operator L the 
eigenvalues problem is formulated as follows 

L | f a ,p)) = W a ,f3 | fa,p)), (7) 

where | f a ,p)) =| ipa >< typ \, w a ,(3 = E a — Ep being the eigenstate and 
eigenvalue of Liouvillian L. The values | ip a > (| ipp >), E a (Ep), are the 
eigenstate and eigenvalue of the total Hamiltonian H, respectively (in this 
case H \ ip a >= E a \ ip a >). The equation ((Tj) can be rewritten by using the 
correlation index v 

L\f v ))=w"\f u )), (8) 

where w u = w a> p and | f u )) =| f a ,p})- in t ne equation ([HD^ = if a = (3 which 
is the set of diagonal operators \ ip a >< ip a \ or "vacuum of correlations" 
and v 7^ is the remaining off-diagonal case. The details of the theory of 
correlations can be found, for example, in the works |14j . [T5] . Eq. (JSJ) is time 
reversible and does not describe the irreversible evolution. Making the well 
known approximations from equation (jSJ) it is possible to obtain the kinetic 
Pauli master equation; however, I will discuss an alternative approach based 
on the new formulation of the eigenvalues problem 

L | = Z] | ((^ | L = ((^ | ZJ, (9) 

where Zj are the complex values, and j is a degeneracy index, since one 
type of correlation index can correspond to different states. In this case, the 
complex eigenvalues Z v - of the Liouville operator cannot be reduced to the 
simple difference of two eigenvalues of the Hamiltonian (as in the case (j7j)) 
and eigenstates | ((^ | are not the multiplications of the wave functions 

(they cannot be obtained from the Schroedinger equation). As it was shown 
by the Brussels-Austin group, this formulation has a broken time symmetry 
and leads to the rigorous description of irreversibility. 

Further development of the approach is based on the theory of subdynamics. 
Subdynamics is a term coined for the construction of a complete set of spectral 
projectors H u 

n" = EI*i»«*il- ( 10 ) 
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The projectors IF satisfy the following relations 

WL = LW ) (commutativity); V'n 1/ = 1, (completeness); 

(11) 

mr =ir<W, (orthogonality). 
Taking fflTT) . it is possible to write down the density operator p in the form 

V V 

where p v = IF/?. Because IFL = ZTF, the operators p 1 ' satisfy the separate 
equations of motion [13] 

i^fi = irur/Kt). (13) 

In other words, dynamical evolution is decomposed into independent " subdy- 
namics"; it is specified by the single parameter v. Operator IF is determined 
by the expression 

rF = (P U + C V )A U (P U + D U ), (14) 

where 

A v = P\P V + D V C V Y X P V (15) 

and C", D u are kinetic operators, P v is a projection operator. Operator C u 
creates correlations other than the v correlations, while D u is destruction 
operator, (operators C u , D u , P u can be found in the works [5], [11], [TJ] and 
in Appendix A). As it has been shown [5], Eq. ( 113"]) can be reduced to the 
Markovian kinetic equation 

^=W(i). (16) 

Here the value i9£ is a non-Hermitian dissipative operator which plays the 
main role in the nonequalibrium dynamics 

d v c = P V L(P V + C V ). (17) 



As it was shown in the work [16], the density operator p can be determined 
by the kinetic operators C u , D u , 

P°p(t) = P° exp{-i$° c t)A°P p{0) + pa C u exp{-i^ c t)A u D u P°p{0). 

The Eq. ([TBI) contains the Markovian part in the first term and non-Markovian 
contribution in the second term. The expression (fTHj) for the density operator 
is essential for our further analysis. 
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In the paper I examine the two-nucleus spin systems which consist of a pair 
of uniform nuclei with a spin of 1/2. The spin systems are determined by the 
H(l) and C(13) nuclei. The interaction between the magnetic moments of the 
nuclei means that the pair of nuclei behave as connected system. For example, 
the transition of one nucleus from the upper level to the lower (the effect of 
the Nuclear Magnetic Resonance (NMR)) can cause the energy transition of 
the adjacent nucleus. In the work the energy transitions are considered to be 
irreversible processes. The two-nucleus spin NMR systems are analyzed in the 
framework of the subdynamics theory. In Section 2 the theory of magnetic 
interaction of the nuclei is discussed. The solution to the complex eigenvalues 
problem is given in Section 3. The time irreversible evolution of the two- nucleus 
spin systems is investigated in Section 4. In Section 5, the results of the work, 
and quantum Zeno effect are discussed. 



2 The magnetic interaction 



Let me examine the pair of the uniform nuclei which are located in the external 
magnetic field, directed along z axis - H z . The Hamiltonian of this system can 
be written down in the form 



Here the first term, Hjh z -> corresponds to the interaction of the system with 
the magnetic field H z and the second, Hnn, is the interaction between the 
magnetic moments of the nuclei. The values I\, I2 are the spin vectors of the 
first and the second nucleus, I z i, I z2 are their z - components, r is the distance 
between the nuclei, n - unit vector in the direction from the first nucleus to 
the second, 7^ is a gyromagnetic ratio 7^ = /i/J, where /x, / are the magnetic 
moment and the spin of the nucleus. 

Further examination of the magnetic interaction of the nuclei will be based on 
the works [T7], [TS]. At the beginning let, the interaction between the nuclei 
be absent. In this case, the projections of the summary spin of the system 
I z = l z \ + I Z 2 on z axis have four different values: 1, 0, 0, -1. The projections 
correspond to the four different states: the spins are directed along the field 
H z (ft); they are directed differently {]/[) or (tl); and the spins are directed 
opposite to the field (O). To describe these four different states, the four 



Hih z +nn — —1nH z (I z i + I z2 
= H IHz + H NN . 



)-7^- 3 [3(/ 1 n)(/ 2 n)-/ 1 / 2 ] 



(19) 
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eigenstates can be determined 



V'3 




(eta), 
(a/3), 



^2 



'04 




((3a), 



(20) 



where the function a corresponds to the spin state (f) (the projection of the 
spin of the single nucleus on z axis is 1/2) and the function corresponds to 
the spin state (|) (the projection of the spin on z axis is —1/2). 
Let me note the validity of the following expressions (where the units h = 1 
and the speed of light c = 1 are used ) : 

1. The action of the spin operator on the functions (i—1 or 2) 



Igi ec 2 ^ ' 



1 i 

Ixi ex. 2 ^ 2 

1 z I 

4i = -at, Iyi P = --a, hi = ~-/3; 

2. The multiplication of the pairs of identical functions 

(aa)(aa) = (a(3)(a(3) = ((3a)((3a) = ((3f3)(f3(3) = 
and different functions 



Furthermore, 



(aa)(a(3) = 0, (aa) (f3a) = and so on. 



I z (aa) = (aa), I z ((3a) = 0(/3a), 
I z (aP)=0(a(3), I Z ((3(3) = -(P(3). 



(21) 



(22) 



(23) 



(24) 



(25) 



Thus, from the expressions fl25l) it follows that the functions (aa), ((3a), (a(3), 
((3(3) correspond to the spin states (tt)> (It); (t-l)> (II); respectively. However, 
only the functions (aa), ((3(3) are the eigenstates of the summary operator 

P = (h+l2? 



I 2 (aa) = 2(aa), I 2 ((3(3) = 2((3(3). 
The functions (a(3), ((3a) are not the same 

I 2 (a(3) = ((a/3) + ((3a)), I 2 ((3a) = ((a(3) + (f3a)). 



(26) 



(27) 



6 



The results ([271) are not satisfactory for the description of the nuclei spin 
system. However, if we substitute the functions ^( (a (3) + ((3a)) and ^((a/3) — 
((3a)) for those of (a(3) and ((3a), then the four new functions are defined as 



0x = (aa), 2 = ((a/3) + ((3a)), ^ = (/?/?), 



v/2 



y/2 



((a(3) - ((3a)) 



(28) 



will be the eigenstates of the operators I z , I 2 with the eigenvalues 1, 0, -1, 
and 2, 2, 2, 0, respectively. From the aforesaid it follows that the pair of 
the noninteracting spins 1/2 in the states 0i, 02, 03 behaves as one particle 
with the spin 1, while in the state 04 as the particle with the spin 0. Since 
the functions 0i, 02, 03, 04 are the eigenstates of the operators I z , I 2 , it is 
convenient use them to describe the interacting spin system. 
Now our aim is to determine the energy of the interacting two-nucleus spin sys- 
tem placed into the external magnetic field H z . The solution to this problem 
can be conveniently carried out in the matrix representation of the Hamilto- 
nian (fT9~j) . Calculating the matrix elements of the operator Hjh z with the use 
of the functions (128]) (in according to the rules (|2T|) - (12^1) ) leads to 



H 



-JnH z 



(\ o o (A 
oooo 

0-10 

\o o o oy 



(29) 



Let me examine the operator H^^. In the Cartesian coordinates the spin 
operators Ii, I 2 and the unit vector n can be represented as 



n = i cos(a) + j cos((3) + k cos(7), 



(30) 



here i, j, k are the unit vectors directed along the coordinate axes and a, (3, 7 
angles between the direction from one nucleus to another and the coordinate 
axes x, y, z. Let me switch over to polar coordinate system. In this case, 

cos(a) = sin(^) cos(yj), cos(/3) = sin(6') sin(<^), cos(7) = cos(6 ) ), (31) 

with 6 and <p being polar and azimuthal angles. Using the expressions (130|) . (13T1) 
we obtain 



h nn =[i zl i z2 - \(i{i+ + A+^h^o - \(ith2 + i z i4h 2 N yi 



(32) 
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where the action "*" is the complex conjugation, 



(I xl ± %l 



yij 



(33) 



and 



Y 
Y 2 



r- 3 {l -3cos 2 (#)), 

r~ 3 sin(6 l ) cos(6 l ) exp(— zy>), 

r~ 3 sin 2 (6 l ) exp(— . 



(34) 



Calculating the matrix elements of the operator H^n leads to the matrix 
representation 



H 



NN 



1 2 



Y 

S^ZY* 
-3Y* 




-3^2^ -3Y" 2 0^ 
-2F 3^^! 
3^1? Y 

o o oy 



(35) 



It has been shown |XT] that the influence of the nondiagonal elements in the 
first approximation can be disregarded since the changes of the energy lev- 
els, caused by them, are small. Taking into account the latter fact and the 
expressions (EES]), fl2U|) . f l3"5"|) for the Hamiltonian H IHz+NN , we have 



H 



IH Z +NN 



-In 



( E Z - \ lN Y Q 

\i N Y Q 












V 



o 
o 



-H z - \ lN Y Q 
OJ 



(36) 



It follows from the expression fl36l) that the interaction of the nuclear spin 
system with the external stationary magnetic field H z as well as interaction 
between the magnetic moments of the nuclei lead to the splitting of the system 
energy levels into the values 



E 1 = - lN H z + -j 2 N Y , E 2 
E 3 = lN H z + - 7 ^y , E 4 = 



■^nYq, 



(37) 



0. 



From the expression (136]) it follows that the singlet state 04 is not affected by 
the interaction since it corresponds to zero summary spin and consequently 
to zero magnetic moment. 

The transition between the energy levels ( 13"T|) can occur as a result of the 
external, variable magnetic field Hi perpendicular to the field H z (the NMR 
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effect) influencing the system. In the simplest case, the components of the field 
Hi change according to the law 

H x — Hi cos(cot), H y — —Hi sin(a;t), (38) 

where u is the frequency of the rotating field in the xy plane. 

In this case, the Hamiltonian of the disturbance Hju xy can be written down 

as 



H iH xy = -j N Hi(I x cos(ujt) - IySm(ojt)), 



(39) 



where I x = I x i + I X 2, I y = I y i + The expressions (139]) and (128!) lead to the 
matrix representation of Hm 



( 



H 



IH X 



' V2 







exp(iut) o\ 



exp(— iut) exp(iwt) 
exp(-iujt) 

o o o oy 



(40) 



Since the transitions are caused by the variable disturbance Hin xy , then it is 
obvious from the expression ([4*01) that the ones are permitted only between 
the states <pi, (p2 and between the states 02, 4>3- m the first case, the transition 
energy Ai? 12 is determined by 



(41) 



A£ 12 = j N (H z + - lN r- 3 (3cos 2 (6) - 1)J, 
in the second case, the transition energy Ai? 2 3 is determined as follows 



AE 23 = lN [H z - - 77V r- 3 (3cos 2 (#) - 1) 



(42) 



The expression (140]) corresponds to the monochromatic model. In order to 
pass to the general case it is necessary to make the following replacements 



Hi exp(iujt) — > J A(tu)Hiexp(iujt)duj, 

— oo 

oo 

Hi exp(— iut) — > J B(iu)Hiexp(—iut)du, 



(43) 



where the coefficients A(u), B(u>) determine the distribution of the radiated 
energies. 

As it has been shown, in the first approximation (136]) the functions (128]) can 
be defined as the eagenstates of the Hamiltonian ( TT9l) . Therefore, I assume 
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that any function which describes the behavior of our spin system can be 
represented as expansion in terms of the latter 

= ai0i exp(-iE x t) + a 2 <\> 2 exp(-iE 2 t) + a 3 3 exp(-iE 3 t) 
+ a404exp(— iE^t), 

where the energies Ex, E 2 , E 3 , E4 satisfy the relationships (I3"T|) . 
For the hermitian conjugated function I have 

0t = a *0t exp^Eit) + a* 2 (plexp(iE 2 t) + exp(i£ 3 t) 
+ al4>\exp(iE 4 t). 

The expressions which were obtained in this section are the basis for the 
further model development. 



3 Complex eigenvalues problem 



Let me examine the transition from the level E 3 (the projection of the spin is 
-1) to the level E 2 (the projection of the spin is 0) as an irreversible process, 
when the energy A.E23 is emitted. The solution to this problem will be carried 
out on the basis of second quantization method. In this case, the functions 
A(oj), B(u) are defined as operators 

A(u) g{u)a ] (u), B{oj) g(u)a(u). (46) 

Here afl(u), a(u) are the operators of creation and destruction of a photon 
with the energy uj. I assume, these operators satisfy the commutation relation 

[tt(u) ) o t ((/)]=%-K/). (47) 

The value g(u) is determined by the fact that its square is the distribution 
function satisfying the condition 

00 

J g 2 (uj)du = 1. (48) 

—00 

The functions (j44p . fT45l) are defined now as operators too. The coefficients 
in the expression f|44|) . where i — 1, 2, 3, 4 specify destruction operators of the 
state with projections of the spin 1,0,-1,0, respectively. For a* in (1431) we must 
make a* — > a\. The values a\ specify creation operators. It is assumed that 
the operators a%, a\ satisfy the commutation relation 

[a m ,a) m ,} = 6mm', m (m) = ±1, 0. (49) 
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The Hamiltonian of the two-nucleus spin NMR system Hnm have to be de- 
termined as 

Hnm = 4 } (Hih z +nn + H IHxy )(j) = H IHz + H NN + H IHxy , (50) 

where the operators in the brackets are determined by the expressions (1191) . (1391) 
(in the expression ( 1501) the previous designations of the Hamiltonians are pre- 
served) . 

In the model, the complete Hamiltonian of the system I define as H = Hq+Hj, 
where H is a free part H = H w + Hjh z +nn and Hj = Hjh xv = XV is a inter- 
action part (with A being a coupling constant). The free photon Hamiltonian 
is determined as follows 

= J ooa) (oo)a(oo)doo . (51) 

— oo 

In this case, \ oo >= oo \ oo >, where the state | oo >= a^(oo) \ vac > is a 
eigenstate of the free Hamiltonian H w and oo is a photon energy ( | vac > is a 
vacuum of the state: a{uj) \ vac >= 0). I determine the states with projection 
of the spin on z axis equal ±1 and 0: | ±1 >= a{ 3 | vac >, | >= a\ A \ vac >, 
respectively. The states | ±1 >, | > are the eigenstates of the Hamiltonian 
Hih z +nn- H IHz+NN I ±1 >= Ei i3 | ±1 > and H IHz+NN | >= E 2 ^ | >. 
Furthermore, another two states | ±1, u >, | 0, oo > correspond to the nucleus - 
photon system. For them, the relationships H \ ±l,oo >= (go+Ex j3 ) \ ±l,oo >, 
H | 0,oo >= (oo + #2,4) \ 0,oo > are valid, where | ±l,oo >= a\^a)(oo) \ vac > 
and I 0,c<; >= a\ tA a)(oo) \ vac >. 

I determine the eigenvalues problem for the complete Hamiltonian H = + 
Hih z +nn + H IHxy for the state that corresponds to the projection of the spin 
on z axis equal to —1. As it has been mentioned above (Introduction: the 
relationships ([1]), (T5])), the description of the transition E3 — > E2 (transition 
I —1 >— >\ 0,oo >), as an irreversible process, has to be executed outside the 
Hilbert space. It is assumed that the eigenvalue Z_i of the Hamiltonian is 
complex. Then I have the equations for the right eigenstate | <£>_i > and for 
the left eigenstate < <J5_i | of the Hamiltonian H 

H I <£>_i >= Z_i I y?_i >, < I H =< <£_i I Z_i. (52) 

According to the method developed in the works [5] and [6], I expand the 
values I <f-i >, < Cp-i |, Z_i in the series 

00 00 
I v . x >= E A" I <pM >, < ip. x |= E A" < ifcl I, (53) 

n=0 n=0 

00 

Z_ 1 = £A n zi n 1 ) , (54) 

n=0 
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where 

><g>=i-i>, <^°j i=<- 1 1 ^ (o) 

A = e — the elementary charge. 



m (V \ — I 1 \ s ,nW |_/ 1 I 7 " _ 
V 5 -! > — I _1 >) < V 5 -! I" < _i I) Z -l - ^3, 



The complex eigenvalue problem (|52|) is solved analogously to the works [6], [T6] . 
In this case, for the coefficients Zl^ it is easy to obtain 

n-l 

=< -l | y | <^T 1} > - E ^-1 < - 1 1 v-T° >, (56) 

i=i 

where the interaction term is determined by the expression 

<-l\V | ^T 1} >= / g{w) exp (-i(^2 - E 3 + u)t) 

eV2 _4 (57) 

x < 0,w | > dcu. 

Using ( JBTl) . for Z_i the following result is valid 

Z-i -E 3 = J g{u) exp {-i{E 2 - E 3 + u,)t) 

x < 0, a; | V9_i > du - (Z_i - S 3 )(< -1 | > -1). 

The determination of the matrix element < 0,u \ (p^\ > follows from the 
obvious relationship 

i OO OO 

< 0,u, |(i^oE en I ^-i > +eVY,e n I ^ >) 

oo n=0 oo n=0 ' m 



n=0 n'=0 

hence it follows 

(n) -1 



< 0, u | > 



i?2 4" w — 

x (60) 

x (< o,w | y | ^T 1} > -E z -i < I v-T° >)• 

The transition | —1 >— y\ 0,oj > is a time ordered irreversible process. There- 
fore, the time ordering of the expression (loUl) has to be introduced. For this 
purpose, it is necessary to make the replacement 

1 1 

E 2 + co-E 3 E 2 + uo- E 3 -ie^ (61) 

Here z/ - correlation index corresponds to the initial state | — 1 > and fi - 
correlation index corresponds to the final state | 0,oj >. In our case, it is 
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possible to introduce the redesignation = e > ( the real infinitesimal 
e^ v determination details can be found in the work [6]). For the first matrix 
element of the expression (IBTij) I have 

< 0, u | V | ^LT 1} >= - J ^9(^) eMKE2 ~E 3 + co)t) 

x < -l i ^t x) > • 

oo 

Substituting the expressions (IBTI) . (|62|) for (|60|) and summing up J] e n result 

n=l 

in 



< 0, w | ip-i >= exp(i(E 2 - E 3 + u)t) 

X < -1 | > — — — -. 

E 2 + lo — Z_x — is 

The expression (1631 leads to 



The formal expression 
1 



(63) 



+ t7c5{E 2 + u-E 3 ) (65) 



E 2 + u — Z-i — is E 2 + uj — E 3 
makes it possible to rewrite the relationship for Z_i in the form 

Z_ x = E 3 - 2 7 -i, (66) 

where E 3 is a renormalized energy of the initial state 

2 tj2 00 i 

2 y 1 J E 2 + u-E 3 (67) 

P being the principal part of the integral 

and 

7-i = ^tfjy (A£ 23 ). (68) 

It is known [19J that the rate of the transition between two states | m > and 
| w! >, where m (m') is the projection of the spin / on the z axis is defined 
by the expression 

W m ^ m > = ^i 2 N Hl(I + m + - m)f(u), (69) 
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where ml = m ± 1, and f(oj) is a distribution function with the width A. The 
distribution function describes, for example, Lorentzian or Gaussian distribu- 
tion of energy. Lorentzian distribution is described by the expression 

= -JvTT 5 = l A ( 7 °) 

71 Z + (U — OJqY 2 



with the norm 



f(u)duj = 1. (71) 



Since the singlet state 04 corresponds to zero magnetic moment and does not 
participate in the energy transitions, we can conditionally assume the spin of 
the two- nucleus system participating in the interaction, / = 1. In our case, 
m——l while m! — 0. Then from the expression fl69|) it follows that 

W-i-o = tt 7 M/(w). (72) 



Comparing the results (J5Sj) . (I72|) we see 

W^-l->.o/a;=cJo = 2 7-lA 2 (^))=/(^o)' ( 73 ) 

where the energy cu = Ai? 23 and (? 2 (wo) = /(^o)- Thus, the expression ob- 
tained for 2 7 _i coincides with the rate (172|) under the resonance condition 
a; = Uq. The conclusion (1731) is in agreement with the results of the work |lbj . 
in which the transition | —1/2 >— >\ 1/2, u > was investigated. It has been 
shown that the expression for 27_x/ 2 coincides with the rate W_ 1/2-4.1/2 under 
the resonance condition ujq = E_y 2 — Ey 2 - in the work [20J, where the irre- 
versible evolution of the unstable n~ - meson was examined, the value 27 p?r 
was defined as the rate of the 7r~ - meson decay. 

Thus, on the basis of the Brussels-Austin group approach, the value of the 
rate is obtained as the solution to the eigenvalues problem, on the basis of the 
complex spectral representation without any mnemonic rules. 



4 Evolution of the two-nucleus spin system 



The description of the irreversible evolution of the two-nucleus spin system 
is carried out on the basis of the subdynamics theory. I examine the matrix 
element 

((-1; -1 I P°p(t))) (the expression ATS]) ): 

((-1;-1 I p(°V(t)>> = (<-!;-! I PW>> m 
=< -1 I pit) I -1 >, 
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where the projection operator P° is determinated in Appendix A. Let me 
represent the operator A (JTB1) as series expansion. In contrast to the work |16j . 
where the opreator A was examined in the first approximation, here I define 
it, taking into account the second term 

A = P°(P° + D°C°)- l P° ^ P° - D°C°. (75) 

In this case, for the matrix element fl74l) it is obvious 



(76) 



(77) 



((-1; -1 | p(t))) = ((-1; -1 | exp(-^)PV(0)}) 

- «-l;-l | exp(-i#° c t)D°C P p(0))) 
+ ]T((-1;-1 \C v e X p(-i^ a t)P v D v P°p(0))). 

For the first term of the expression (176]) I have 

«-l;-l | exp(-i^t)P°p(0))) 

= -1 | exp(-i^t) | -1; 

where the initial density operators was taken as 

p(0) =| -1 >< -1 |, i.e., ((-1; -1 | p(0)» = 1. (78) 

Using the determination of the dissipative operator ^ (Eq. ffTT|) ) and the re- 
sults of the work [TB] , in which the calculating details the analogous expression 
for a one-nucleus spin system are given, it is easy to obtain 

((-1; -1 | exp(-z^)P°p(0))) = exp(-2 7 _ 1 t). (79) 

Now let me examine the second and the third terms of the expression (176 p . 
For the second term I have the result 

((-1;-1 | exp(-^)P°C P p(0))) 

. n ^l 2 N Hl( ? f{cu) \ (80) 

= exp -2 7 _it / JK ! , - duj + c.c. . v ; 

2 \J [CO — u — IE) 1 J 

The expression (|80|) obtaining details are represented in Appendix A of this 
work. Further, for the third term, the following expression is correct 

£«-L;-l | ^exp(-^)P^P°p(0))) 

1n h I ( 1 ff x exp(-i(w - uo)t) 

f(u)—, —^-du + c.c. 
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The integration in the expressions (180)) . ( IHTj) is carried out in the investigated 
kinematic region, with the function f(oo) satisfing the condition 



f(u)du; = 1. (82) 



Comparing the results (I80I) and (IHTj) at the moment t = 0, one can see 
«-l;-l | exp(-<t) J D C P p(0))) 4=0 

= -1 | ^exp(-^)P^P p(0))) t=0 ( 83 ) 

whence the correct value for the matrix element at the initial moment of the 
time follows: < -1 | p(0) | -1 >= 1. 

Taking into account the relationships (j7ij) . ([76]) and (j?5j). ([50]) . (jglj) I obtain 
the expression 

< -1 | p(t) | -1 > 

= exp(-2 7 -it) - exp(-2 7 -it)4^-f / 7 ^ : + c.c. 

2 V (w-wo-^) 2 / (84) 

, l%Hl ( J ft ^ exp{-i{u - u )t) 
2 w (a; — w — *e) 

The integrate calculation features the expression (18^1) are examined in Ap- 
pendix B. Calculating the integrals finally results in 

< -1 | p(t) | -1 > 

= exp(-27_it) + exp(-27_it) 2 j^ 2 ^ {25 + ttu + 2u arctan(y )J 

+ 7^8 W) sin ^) + ^(t) cos(w *)) - exp(-5t)^^. 
OITUqI o z 

In the expression (18 5 j) 

^)^/\ex P (- ^V( 1 ^)-g( 1)+ 66)f + 6a i (86) 



/• , -a 7 (2 + 46)£ + 4a 5 6£ 3 

= y deexp(-o — ^ — ^ -. (87) 

D(0 = a 8 (l + 6) 2 + a 6 (2 + 6(2 + 46))£ 2 + a 4 (l + b(-2 + 66))£ 4 
+ a 2 6 (_ 2 + 46)e 6 + 6 2 e 8 ; 1 J 

a = and b = (uq/5) 2 . 
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5 Results and discussion 



In Fig. 1, 2 the transition from the energy level E 3 to the level E 2 of the pair 
of the connected hydrogen nuclei H(l) (mass number A = 1) is investigated. 
The transition is accompanied by an uq = AE23 energy quantum emitting. 
The evolution of the density matrix is examined for a long time scale (Fig. 1) 
and for a short time scale (Fig. 2). It is shown that the behavior of the com- 
plete density matrix f[8"5]) (solid line) is slower than the exponential one (1791 
(dashed line). In the asymptotic time region, t — > +00, it approaches to the 
latter. The complete density matrix is determined on the short time scale, 
where the deviation from the exponential, Markov evolution is seen (Fig 2). 
The behavior of the nuclear spin system is found depending on the initial dy- 
namic parameters selection. For example, Fig. 1 presents the calculations for 
the field Hi = 37 Oe. In this case, the deviation from the exponential evolution 
is not considerable. Fig. 3,4 present the calculation of the time evolution of 
the connected carbon nuclei C(13) (mass number A = 13). The calculation is 
carried out for the field Hi = 100 Oe and for the field Hi = 150 Oe (Fig. 3). In 
this case, the behavior dynamics of the nuclear spin system coincides with the 
foregoing case. The interaction with the field Hi = 100 Oe leads to a noticeable 
deviation from the Markov evolution. In the case Hi = 150 Oe the deviation 
is not considerable. On the short time scale, (Fig.4) the nuclear spin system 
behavior is slower than in the H(l) case. In both cases, the calculations were 
made for: 5 = 10 6 sec -1 , 9 = 30°, r = 2 x 10~ 9 cm. The magnetic moments of 
the nuclei have the values: 2.7927 for H(l) and 0.702381 for C(13) (the values 
of the magnetic moments are given in the units of the nuclear magneton). 
The suppression of transitions between quantum states obtained on a short 
time scale in the course of the present work may be thought of as the so- 
called Zeno effect occuring as a measurement result in quantum physics. At 
present, Zeno effect is discussed widely in the literature (see, for example, the 
works [23]- [22] )• Sudarshan and Misra's pioneer work [22] showed that con- 
tinuous observation of the decay process makes the latter impossible as the 
observation influences the decay. In other words, the decay slows down under 
the observation and in the case of continuous, ideal measurement it becomes 
impossible. Thus, Zeno effect contains a subjective, anthropomorphic element 
as quantum theory describes the fact that we can speak about a quantum 
world after taking measurement rather than quantum world itself. The result 
depends on the observer and on observation conditions. This paradoxical re- 
sult was subsequently examined in connection with the problem of quantum 
transitions in atoms [23] . The theoretical aspects of the paradox are discussed 
in the monograph [20] (see also Khalfin [27] and the work [28]). 
Quantum Zeno effect was investigated by the Brussels-Austin group in the 
works [29] - [51] . In the work [29], the system of harmonic oscillators coupled 
with an external field was examined. Zeno effect was discussed in connection 
with the problem of quantum decoherence. It has been shown that quantum 
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Zeno time serves as a lower bound for the decoherent process. In the work |30j . 
a particle coupled with a field in the thermodynamic limit on the basis of the 
complex spectral representation of the Liouvillian was discussed. The results 
show the time derivative of the photon number density asymptotically van- 
ished at t — > 0. It was shown that the influence of Zeno effect is noticeable and 
determines the subsequent system dynamics. The short-time behaviour of the 
survival probability in the frame of the Friedrichs model was investigated in 
the work [31]. The problem was associated with a photodetachment process, 
quantum dot, and a hydrogen atom. It was discovered that Zeno effect occurs 
in a very short time region. 

In connection with the quantum computation problems, the studies dealing 
with creating quantum systems for which Zeno time region is long and which 
make quantum calculation possible are of great importance. In this connec- 
tion, the work of Center for Quantum Computation, University of Oxford 
is noteworthy. The work [22] presented the preliminary experimental results 
demonstrating the quantum Zeno effect in NMR for H(l), C(13) nuclear spin 
systems. The specific features of the future NMR experiments, in which Zeno 
effect can be demonstrated, are discussed. The authors state, that " Zeno effect 
in NMR experiments has not previously been explored, probably as a result 
of the dificulty of performing true quantim measurement." 



6 Conclusion 

The investigation of the evolution of the two-nucleus spin systems interacting 
with the magnetic field in the frame of the complex spectral representation of 
the Liouvillian gives the rigorous description of the irreversibility. The short 
time (Zeno) region determines the memory effect on the quantum level. In 
the approach, Zeno effect is considered as an objective dynamic process tak- 
ing place in a quantum systems on a short time scale without using anthro- 
pomorphic elements. The study of the two-nucleus spin systems reveals the 
significant role of the effect which, however, depends on the characteristics of 
the magnetic field interacting with the nuclei. The investigation of the relax- 
ation process for the two-nucleus spin systems demonstrates the utilization of 
the Brussels-Austin formalism in a practical problem that may be tested in 
experiment. 
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Appendix A. The second approximation 

Let me examine the second term of the expression (1761) 
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((-1; -1 | exp(-i$° c t)D°C o P o p(0))). The operators C°, D° have the form 

D° = y P°Lj P». (A.2) 

In our model, the projection operators P°, P 11 (fi 7^ 0) are determined as 
follows 

1 °° 

P° =| -1; -1»«-1; -1 I +- J I u;0w))((0 w;0w| du, (A.3) 



OO 



P^= / I -l;0w))((-l;0w I du + / I 0w;-l))((0o;;-l | du, (A.4) 



—00 



where the value f2 is designated so that ^ = 1 and 5(0) is the delta-function. 
The expressin ( ISUj) is obtained in the limit Q — > 00. In this case 



((-1;-1 I exp{-iti° c t)D C P p{0))) 

+ i£ (A. 5) 



"I I exp(-z^) I -1; -1 I Li 1 P" 



x-pVL/|-l;-l»«-l;-l|p(0)». 

For the first cofactor of (1A.5j) . I have the result (1791) . For the second cofactor, 
taking into account the determinations 

Lj = Ht x 1 - 1 x Hi, (Ax B)C = ACB, (A.6) 

((a; /3\Li\ a'; /?')) =< a \ Hi \ a' > Spp - ^ < P' I I >, (A.7) 

where 8pp (5 aa /) is delta-function if the indices (3', (3 (a, a') correspond to the 
continuous spectrum and it is Kronecker's symbol if the indices correspond to 
the discrete spectrum, it is possible to write down 

— + ie — 

_ vV< _1 I #1 I P >< P I Hi I _1 > 
= V V (P p - P 3 - 

< p I if/ I -1 >< -1 I P/ I p > 
+ (E 3 -E p -iey 

7^? / 7 /(*>) 



(A.8) 



(oj — LU — is) 2 



duj + c.c. 
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Here, index p corresponds to the state | 0, uj >, E p = E 2 + oj, uiq = E 3 — E 2 . 
In (1A.8[) . the sum determines the summation (integration) over all discrete 
(continuous) indices. The results ( I75|) and ( ]A.8j) lead to the expression 

Appendix B. Calculation of the integrals 

I examine the second integral of the expression (jH 



„, . exp(— i(co — u)o)i) , 

f(u)-p^ r^du. B.l 

(to — ujq — ley 

Since, t > 0, we have to deform the contour of integration Tr in the lower half 
plane [22]. The application of the residues theorem leads to 

J exp(— i(z — Lo )t)F(z)dz = —2iri ^ Res [exp(— i(z — u) )t)F(z)], 2) 



where 

F « - v^k? (R3) 

and the sum comes from the poles zi of the function F(z). In the expres- 
sion flR2|) 



R 

exp(— i(z — Uo)t)F(z)dz = J exp(—i(co — u )t)F(co)dco 

° o (B.4) 
+ J exp(—i(z — uo)t)F(z)dz + J exp(—i(z — uo)t)F(z)dz. 

1r -iR 

If f(z) is the Lorentzian distribution ([70]) then, 

lim / exp(-izt)F(z)dz -> 0. (B.5) 

IR 

In the lower half plane the function F(z) has one pole z = ojq — id. In this case 
-2m £ Res[exp(-i(z - u )t)F(z)\ = _ ex P(~^) , (R6) 

Z = Z; 

Using the result obtained in the limit — >■ 00 we can write down 



00 —too 







exp(— i(ou — u )t)F(u)duj = j exp(—i(z — ujo)t)F(z)dz — -. (B.7) 
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The replacement z = —iy, makes it possible to reduce the result (1B.7j) to the 
form 



exp(— i(u — uo)t)F(u)du = —i exp(—i(—iy — uo)t)F(—iy)dy 

o ' (B.8) 
exp(— 5t) 

I 2 ' 

Then the replacement y = £/t, where £ > leads to the final expression 

^-^-^Y (b.9) 
= -=-^-(A(t) sin(w t) + B(t) cos(w„t)) - 2 exp [~^ , 

It is necessary to note, that the result for the integral (IB.lj) in the foregoing 
work [16] was obtained without taking into account the residue ( 1B.6I) . This 
was made because the role of the latter in the investigated kinematic region 
was not essential. 
The integral 

7 ^ . r^ du (B.IO) 

[U — U — IE) 1 ' 
u 

can be obtained from the expression (IB. 71) for the case t = 

oo — ioo 

jF{u)du= J F{z)dz-—. (B.ll) 



Using the "Maxima" program and taking into account the complex conjugion, 
I then have the result 

7 /M . j_ 25 + 7ru; + 2a;oarctan(f ) 

du + c.c. = ^ — . (B.12) 



(uj — loq — ie) 2 ti5 2 ujq 
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Figure Legends 



Fig.l: H(l), H z = 10 4 Oe, H x = 25 Oe and H x = 37 Oe; solid line is for 
complete density matrices f[55"j) . dashed line is exponential one ( 175|) . 

Fig.2: H(l), = 10 4 Oe, Hi = 1 Oe; the designation of the lines as in Fig. 
1. 

Fig.3: C(13), H z = 10 4 Oe, H x = 100 Oe and H x = 150 Oe; the designation 
of the lines as in Fig. 1. 

Fig. 4: C(13), H z = 10 Oe, Hi — 1 Oe; the designation of the lines as in Fig. 
1. 
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